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ABSTRACT 

Sections  1-6  give  the  solution  of  the  main  problem  for  a  spheroidal 
earth  with  the  potential  limited  to  the  principal  term  and  the  second  har¬ 
monic  which  contains  the  small  factor  kg.  The  solution  is  developed  in 
powers  of  kg  in  canonical  variables  by  a  method  which  is  basically  the 
same  as  that  used  in  a  different  problem  by  von  Zeipel  (1916).  The  per¬ 
iodic  terms  are  divided  in  two  classes:  The  short-period  terms  contain 
the  mean  anomaly  in  their  arguments;  the  arguments  of  the  long -period 
terms  are  multiples  of  the  mean  argument  of  the  perigee. 

The  periodic  terms,  both  of  short  and  long  period,  are  developed  to 

o 

O(kg);  the  secular  motions  are  obtained  to  0(kg  ).  The  results  are  ob¬ 
tained  in  closed  form;  no  series  developments  in  eccentricity  or  inclina¬ 
tion  arise.  The  solution  does  not  apply  to  orbits  near  the  critical  inclina¬ 
tion,  63.4°,  but  is  otherwise  valid  for  any  eccentricity  <  1  and  any 
inclination. 

Section  7  gives  the  long-period  terms  and  the  additions  to  the  secu¬ 
lar  motions  caused  by  the  fourth  harmonic  in  the  potential;  section  8  gives 
the  contributions  by  the  third  and  fifth  harmonics;  section  9  contains  for¬ 
mulas  for  computation. 


1.  The  equations  of  motion.  The  equations  of  motion  of  a  small  mass 
attracted  by  a  spheroid  are 
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A  9U  d2y  _  8U 
dt2  9x  dt2  By 


d2z  _  8U 
dt2  8z 


with 


(1  -  3  sin2)3) 


o 

+  — -  (1  -  10  sin  |3 


sin^/3) 


The  equatorial  plane  of  the  spheroid  is  taken  as  the  xy  plane;  0  is  the  lati¬ 
tude,  and  if  M  is  the  mass  of  the  spheroid  and  k  the  Gaussian  constant, 

M  =  k2M. 

The  main  problem  of  artificial  satellite  theory  may  be  considered 
the  problem  with  k^  =  0.  Let  I  be  the  instantaneous  inclination  of  the 
orbital  plane  with  the  equatorial  plane,  g  the  distance  of  the  pericenter 
from  the  ascending  node,  and  f  the  true  anomaly.  Then 


sin  j3  =  sin  I  sin  (g  +  f)  , 

2  sin2/3  =  sin2  I  [  1  -  cos  (2g  +  2f)  ]  , 
and  the  disturbing  function  can  be  written 

R  =  ^2  /  1  +  3  cos2j\  +  /3  _  3  cos2j\  a^cos(2g+2f) 

a3  L  \  2  2  J  r 3  \2  2  )  r3 


Let  a,e  be  the  osculating  semi -major  axis  and  eccentricity,  respec¬ 
tively.  The  Delaunay  variables  are  then 

i  =  mean  anomaly 
g  =  argument  of  the  pericenter 
h  =  longitude  of  ascending  node. 


L  =  (Ma)1/2 
G  =  L  (1  -  e2)1/2 
H  =  G  cos  I 
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With  these  variables  the  equations  become 


dL 

9F 

di 

_  9F 

dt 

TJ 

dt 

dL 

dG 

8F 

dg 

9F 

dt 

~  9g 

dt 

9G 

dH 

ii 

1  07 

I  *3 

dh 

_  9F 

dt 

9h 

dt 

9H 

3HVCOs(2g+2f) 
2  G  V  r3 


(2) 


The  parts  of  F  not  exhibited  in  terms  of  the  Delaunay  variables  may 
be  expanded  in  Fourier  series  as  follows 


—  cos  (2g  +  2f)  =  £  Qj  cos  (2g  +  jl)  *  a2 
r3  j=  “°° 


The  coefficients  P.,  Q.  are  power  series  in  the  eccentricity  e;  the  lowest 
power  of  e  is  e^  in  P^;  e 1 ^’2I  in  Qj,  j  /  0.  Since  e2  =  1  -  G2/L2,  the  de¬ 
rivatives  of  a  function  \p(e)  with  respect  to  L  or  G  may  be  obtained  by 


9^  - 1 
8L  e  9e  ^3 


dip 

9G 


1  G 
e  9e  l2 


) 


(3) 
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The  functions  p  that  will  arise  are  functions  of  a/r  and  f.  Well-known 
formulas  of  elliptic  motion  are 


8e 


cos  f, 


(4) 


8f 

9e 


sin  f . 


(5) 


Some  important  properties  of  the  Hamiltonian  F  may  be  noted: 


(a)  the  independent  variable  t  is  not  explicitly  present  in  F,  hence 
the  integral  F  =  constant  exists; 

(b)  the  variable  h  is  not  present  in  F; 

(c)  the  coefficient  Qq  in  is  zero. 

The  value  of  this  coefficient  is  the  constant  term  of  the  expansion  in 
terms  of  the  mean  anomaly  of  a^r’^  cos  2f:  hence 


%-f 

*  0 


it  a3 


cos  2f  dl . 


By  making  use  of  the  integral  of  areas  in  the  form 


di  =  k  £_2  <jf 
Ga2 


this  may  be  written 

_  lLr^a  1  T  3  71 

%  =  -  J  -  cos  2f  df  =  _  _  (  (1  +  e  cos  f)  cos  2f  df  =  0  . 

11  G  0  r  g3  o 

This  derivation  may  be  generalized  to  show  that  the  constant  term  of 
(a^/r^)  cos  qf  (p  and  q  positive  integers),  expanded  in  terms  of  the  mean 
anomaly  is  zero  if  q  >  p  -  2  ^  0. 
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2.  Outline  of  the  method  of  solution.  Consider  a  Uansformation  from 
the  variables,  L,  G,  H,  A  ,  g,  h  to  new  variables  L’,  G’,  H’,  J[\  g’,  h\ 
with  the  aid  of  a  determining  function  S  (L',  Gf,  H',  l,  g,  h).  Then,  if 


L 


as  P_as 

aJt  9g 


H 


M.jl- 

as 

,  as  Uf 

,  g  =  —  .  h' 

as_ 

(6) 

8h 

8L* 

9G' 

8Hf  » 

the  equations  in  the  new  variables  will  be 


dL’ 

9F* 

ii 

T3 

8F+ 

dt 

aiT 

dt 

9Lf 

dG* 

9F* 

dg;  = 

9F* 

dt 

dt 

aGf 

dHf 

9F* 

dhf 

9F* 

dt 

9h? 

dt 

9Hf 

with 

F*  (L’,  G',  H’,A  ’,  g',  -)  =  F  (L,  G,  H,Jt,  g,  -)  (8) 


The  dashes  in  the  places  for  h?  and  h  are  used  to  indicate  the  absence  of 
these  variables. 

The  problem  would  be  completely  solved  if  a  determining  function 
were  found  such  that  F*  is  a  function  of  the  variables  L’,  G’,  Hf  only. 

The  differential  equations  (7)  show  that  then  Lf,  G’,  H*  are  constants, 
while  A’,  gf,  h*  are  linear  functions  of  time.  Substitution  of  this  solution 
of  the  primed  quantities  into  (6)  then  yields  expressions  from  which  the 
original  variables  may  be  obtained  in  terms  of  L',  G?,  H’,  X\  g’,  h'  and 
therefore  in  terms  of  t  and  the  constants  of  integration. 

In  the  present  problem  it  is  more  convenient  to  choose  the  determin¬ 
ing  function  S  in  such  a  manner  that  X '  is  not  present  in  F*,  while  g*  is 
permitted  to  appear.  If  this  is  accomplished,  L*  and  H*  will  be  constants, 
and  the  system  is  essentially  reduced  to  one  of  one  degree  of  freedom: 


5 


dG*  _  £Z*  -  .  iZ! 

dt  9g?  dt  9G’ 

in  which  L*  and  H’  are  present  as  constants.  After  this  system  is  solved, 
JL',  h’  are  obtained  by  quadratures  from 

dJt'  _  9F»  dh'  _  9F» 

dt  9L’  ’  dt  9H* 

Preferably,  a  canonical  transformation  by  the  choice  of  a  suitable  deter¬ 
mining  function  leads  to  an  equivalent  result. 

3.  Solution  to  the  first  order.  The  determining  function  is  obtained  by  a 
method  used  by  H.  v.  Zeipel  (1916)  in  a  qualitative  study  of  the  motions 
of  minor  planets.  In  this  procedure  a  development  in  powers  of  the  small 
parameter  kg  is  introduced.  If  Fj  is  written  for  R,  the  original  Hamil¬ 
tonian  is 


in  which  the  subscript  denotes  the  power  of  kg  present  as  a  factor.  Fq  is 
a  function  of  L  only.  Let  also 

s  =  s0  +  Sj  +  S2  . 

F *  =  F0*  +  Fj*  +  F2* . 

It  will  be  convenient  to  choose 
S0  =  Vi  +  G'g  +  H'h. 

If  then  the  expressions  for  L,  G,  H,i  g',  h',  given  by  (6)  are  sub- 
stituted  into 


F  (L,  G,  H,i  ,  g,  -)  =  F*  (L\  G',  H',  g',  -)  , 


6 


there  results 


^  8i  9g  8h  J 


-V*  V  j|.  • 


+  F.*  (l1,  G\H',  ii  , 

1  \  8G' 


Expanding  everywhere  by  Taylor’s  theorem  to  the  second  power  of  l^j 
the  result  is 


F0  (L1)  + 


9F0  es,_  3F0  3S2  J  32Fq  /SSj-  2 


+  — 


8L'  3 1  3L'  9 JL  2  8L>2 


3F-.  3S,  3F,  3S, 

+  F,  (L',  G',  H',  j? ,  g,  — )  +  — -  -f  + - 

1  8L’  8i  3G'  3g 


=  F0*  +  Fj*  (L'.G’.H',  -,g,  -)  + 


3Fj*  8SX 


3g  8G' 


+  F»*  (L',  G',  H',  -,g,  -)• 


Parts  of  corresponding  order  in  k,  on  both  sides  yield: 


order  0: 


F0(L')  =  F0*(L') 


(9) 


order  1: 


order  2: 


3F0  9Sj 


+  F,  =  F,  * 


8L-  91  "  ‘‘  '  1 

8F0  8S2  1  82F0  /3S1^2  +  ^1^1 

+  8L’ 

+  8Fj  3Si  _  F  t  +  8Fi* *  »Si  _ 


(10) 


8L'  3l  2  0L,2  \3i 


3G'  3g 


3g  8G' 


(11) 
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in„  could  be  carried  on  Indefinitely,  but  for  current  practical 
requirements  consideration  of  parts  beyond  the  second  order  in  k,  ap¬ 
pears  to  be  unnecessary.  Put 


*1  =  Fls  +  F1P  ’ 


•  is  the  part  independent  of  i ,  Flp  the  part  dependent  on  i  • 

in  which  Flg  is  the  pan  ip  ^  ^  gecular  part  o{ 

in  the  terminology  ^^  Zdic^rt.  Hence,  if  for.  the  sake  of 
the  disturbing  function,  Flp  the  period  f 

brevity  the  notation 


A  =  -  -  +  —  — ~ 

2  2  q2 


3  3  H 

B  =  +  — - — 

2  2  q2 


(12) 


is  used, 
F 


p4k2 


18  l3g3 


A,  F 


/k2 


lp 


(Acrj  +  Bcr2) 


Now  equation  (10)  is  split  up  into  two  equations 

y 

9  jC  T  t  3 


„4k 

8S,  ^2  /*  ro  1  F  *  -  -  —  A 

1-  Z  (Aa1+Ba2),  Fx  - 


(13) 


in 


which  use  is  made  of 


9L’ 


By  integration  may 


be  obtained  in 


the  form  of  an  infinite  sen  a; 


8 


Si  = 


p\ 

,.3 


00  +°° 

A  £  -A  P.  sin  ji  +  B  £  A  Q.  sin  (2g  +  ji) 
1  j  3  -oo  j  J 


.  (14) 


No  constants  of  integration  are  required,  since  only  partial  derivatives 
of  Sj  with  respect  to  Lf,  Gf,  Hf  i ,  g,  h  will  be  needed. 

For  the  following  developments  it  will  be  advantageous  to  obtain  a 
closed  expression  for  S^.  The  method  used  at  the  end  of  Section  2  for 
obtaining  the  integral  of  functions  of  the  type  a^r*15  cos  qf  can  be  used  to 
obtain 


cigdJ  =  /  —  cos  (2g  +  2f)  d 4 

3 

r 


G  L 


A  sin  (2g  +  2f)  +  —  sin  (2g+  f)  +  -  sin  (2g+  3f) 
2  2  6 


Hence 

S 


1 


P2k2  f 

=  - £  J  A  (f  + 


m3 


+  B 


(f  -i  +  e  sin  f) 


1 


(15) 


sin  (2g  +  2f)  +  -  sin  (2g+  f)  +  —  sin  (2g+  3f) 
2  2  6 


>  > 


as, 


M2fc2 


9g  rr  3 


B 


cos  (2g  +  2f)  +  e  cos  (2g  +  f)  +  _  cos  (2g+  3f) 

3 


.(16) 


Introduce 


*2 


2. 

M  k2 
,-4 


9 


L  =  L'  + 


=  LW1  + 


_  1_  +  3  IT 


a3  L'3 


2  2  Qi2  /  \r3  G'3 


+  (i  .  1  SLW  cos  (2g+  2{) 


2  2  q,  2J 


G  =  G'  + 


=  GW1  +  y 


L'4  / 3  3  H2 


,,4  \2  2  g,2 


cos  (2g  +  2f) 


+  e  cos  (2g  +  f)  +  —  cos  (2g  +  3f) 
3 


H  =  H\ 


Sj  was  obtained  as  a  function  of  G',  H,  e,  f,  g.  Hence  L'  is  present 
through  e  and  f  only.  To  obtain  aSj/aL’  it  is  convenient  to  obtain  first 


as.  f  r  af 

-  =  - J  A  (1  +  e  cos  f)  —  +  sin  f 


+  B  cos  (2g  +  2f)  (1  +  e  cos  f) 


+  —  sin  (2g  +  f )  +  i.  sin  (2g  +  3f) 
2  6 
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in  which  use  was  made  of 


cos(2g+2f)  +  lcos(2g+f)+  JLcos(2g+3f)  =  cos  (2g+ 2f)(l+ e  cost). 
2  2 


With 


8f 


(1  +  e  cos  f)  — 
9e 


a_  +  JL  ]  sin  f 
r2  L2 


there  results 


as 


1  f*  k2 


h. 


9e 


,3 


A  I  —  —  +  -  +  l)  sin  f 


>r2  L’2  r 


+  -B 


L  Jl,  Oil  -  —  +  l]  sin  (2g+  f) 


2  T  ,2  r 
r  L 


a2  G'  .  a  x  M  sin  (2g  +  3f) 


+  l  r2  L’2+‘r  +  ' 


With  the  aid  of  the  relations  (3)  there  follows  at  once 


1=  X’ 


as 


9L’ 


=  r  -  ii  xl 


^2  LT 


!  3  jT 


+±+  Asin« 


2  2  q,2; 


2  t  1 2 

r  L 


3  H‘ 


4  4  Qt  2, 


L  ^  9L.  -  —  +  l\  sin  (2g+  f) 


r2  L’2 


1 


a2  G'2  a  , 

+  i  —  —  +  —  +- 

r2  L'2  r  3 


sin  (2g  +  3f) 


(20) 
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g  =  g’  ' 


aSj 

aG' 


.  y2  L'2 
g  +  — 


e  G'2 


!_  J.  +  _3  h; 

2  2  g,2 


a  9L  f  —  +  1]  sin  f 
r2  L’2  r 


—  —  -  1  +  1^  sin(2g+  f) 


+  f—  sin  (2g  +  3f) 

Vr2  L'2  r  V 


y2 


,4 

i 

vi 


£  2  +  2L_1  (f  -X  +  e  sin  f) 


2  2  g»  2> 


9  _  is  h; 

2  2  q*  2, 


1  sin  (2g+  2f)+  -sin(2g+  f) 

2  2 


+  z.  sin(2g  +  3f) 
6 


(21) 


h  =  h'  - 


aSj 

aH 


=  h*  -  3r2 


,4 


.,4 


f  - 1  +  e  sin  f 


-  I  sin(2g+  2f)  sin(2g  +  f)  -  -  sin(2g+  3f) 
2  2  6 


(22) 
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In  the  calculation  of  the  coordinates,  A  and  g  are  not  needed  separately 
but  only  the  sum: 


,  «,  .  y2  fh'2  L 

JL+  g  =  x  +  g'  +  —  ' 


e  VG'2  G’ 


1  3  H‘  \[a.2  G’2  a 


-  -  +  — 


2  2  g,2J\t2  l,2 


+  — +  1|  sinf 
r 


+,1-1* 


4  4g,2 


—  +  1  ]  sin (2g  +  f) 

L\r2L’2  r 


+  (—  +  —  +  i.)  sin  (2g  +  3f) 

.r2  L'2  r  3 


y2 


,4 

i 

^4 


_  3  +  is  h; 

2  2  q,2 


({  -JL  +  e  sin  f) 


+  ,  i-HE 

2  2  g,2 


1  sin(2g+  2f) 

2 


+  1  sin(2g+  f)+  —  sin (2g+  3f) 
2  6 


(23) 


Since  (1  -  e2)"1  -  (1  -  e2)  ^  is  divisible  by  e2,  the  first  part  of  the  dif¬ 
ference  (1  +  g)  -  (1*  +  g’)  has  e  as  a  factor  and  not  as  a  divisor. 


4.  Second-order  contributions  to  the  Hamiltonian.  For  use  in  the  evalu¬ 


ation  of  F2*  it  will  be  convenient  to  put 


,3 


p2 


.,3 


cos(2g+  2f)  +  ecos(2g+  f)  +  cos(2g-i-  3f) 

3 


(24) 


so  that 


9S1 

—  =  - —  Bp  . 

3  2 


8g 


L" 
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Also,  define  ,  Tg  by 


1  Q  r«4  T  ^ 

t,  =  1  — 1  =  3  cos  f  -  3  iL  , 
1  e  9e  e  4  q5 


1  9<j2 


1 


t2  =  ~ 

^  e  9e  e 


1  a 


a3  L2 


2  r4  r3  G2 


cos  (2g  +  f) 


1 

+  - 


5  a_ 

2  4 
r 


a3  L2 
+  r3  G2 


cos  (2g  +  3f), 


which  are  easily  obtained  with  the  aid  of  the  expressions  (4)  and  (5). 
Then 


9F 


Ip  = 


91/ 


n\ 

x.7 


G'2 

-  6  (Ao^  +  Ba2)  +  (Atj  +  BTg) 
L'2 


(26) 


8F 


IP  -  . 


8G' 


4. 

P  k2 

7 7T 


*  (Ar1+Br2)+3il'  ll^) 
L1  1  1  G’  n,2  1  1 


(27) 


If  the  terms  dependent  on  JL  are  not  needed  to  the  second  order  in  kg,  the 
only  part  of  interest  in  (11)  will  be 


1 

2 


92F, 


9L 


,2 


91/ 


dSx  8FX 
9  i  +  9G’ 


) 


where  the  subscript  s  designates  the  parts  independent  of  i.  The  contri¬ 
bution  to  the  right-hand  member  that  has  aF^*/9?  as  a  factor  vanishes 
because  F^*  is  a  function  of  LT,  G’,  HT  only. 

The  following  products  are  easily  obtained: 
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1  8  V0S1 


2  3L' 

,2 \Bl 

) 

8Flp 

asx 

9L’ 

8i 

L'10 

!!ib 

3SX 

P8k22 

8G1 

3g 

L,l° 

8Fls 

8Sj 

CM 

CM 

8G* 

9g 

L,l° 

6,,  2 

3  ^  k2 

2 


(Aal  +  b<72) 


-6  (Aaj  +  Bcr g  i  2  +  —  (Aaj  +Ba2)  (Ar1  +  Bt2) 

L'2 


,2 


,L’  H 

G’  g,2 


B  (-ajP2  +  a2p2)  -2L  (ABtiP2  +  B2^) 


T  A  f  ij2 

-3±i-  (AB  +  B— )P 

,4  2 


and  hence, 


6.  2 

*  k2 

1  "  Too 


„2 


+  AB 


2  L’2 

0 
G 


*9»i®2 + 


,2 


(V2  +  Vl) 


r*  *  T 

-  tiP2  -3  —  P2 


+  B2  [-\°i  *—2°l72-^riP2 

L 


+  B 


H‘ 


m2 


^  (-3ctiP2  +  3<72P2)  ‘3  ^P2 


,4 


,4 


s  .(28) 
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The  parts  independent  of  1  may  be  found  by  evaluating  integrals  of  the 
same  type  as  those  which  occurred  in  the  preceding  section.  The  only 
exceptions  are  the  integrals 

JL  f  cos  f  di  =  -e^ 

7T  0 


i  f1  cos  2f  dJL  =  1 
v  <T  .2 


r3  r2 
2  —  -3  _  +  1 


7T  0 


cos  3f  dJL  =  - 


r*3 

2  G_  -  3  _  +  1 


+  3c1 


The  separate  parts  of  F2*  are  as  follows: 


constant  part: 


9  „  2 

“  al 

2  1 


—  —  +  ^  +  —  — 


16  g5 


2  g6 


8  q7 


315  1/ 
16  G9 


a  T  _  +  15  L.5  .  3  L6  .  H5  L.7  +  315  L9 
+  L2  16  G5  G6  8  G7  16 


Gl 


Sum  =  -A^+  + 


4  G5  2  G6 


4  G7 


la2=-9!  +  135  315  L9 

2  2  32  ^5  16  g7  32  r-.9 


G' 


G 

+  —  a2T2 
L2 


+  15  L 5  .  105  L7  +  315  L9 
+  32  g5  16  g7  32  G9 
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G  .  2  L5  .  5  L_7 

'  I  T2p2  =  +  3  g5  2  g7 

7  L5  5  L7 


L  L5 

+  3  -  °2P2  '  ‘  5  +  2  r7 

G  G  G 


coefficient  of  cos  2g: 


9  L  io  L 

-9(7,0,,  =  -  -  —  +  18 

12  4  g5  g7 


63  1/ 
4  G9 


L7  63  L9 


r2  ,  ^  &  L  -  14  _  +  £2  ii- 

+  -<V 2+Vl>  f7T5  ,7  4  rt 


G  5  L5  21  L7 

*T  T1P2  "  +  7  _5  '177 

Li  G  G 

7 

i 

/7 


-3^p2  = 


i  L' 
+  1  — 


G 


^2-  -r 
L  +  G 

G 


-2-£ 

G 


L 

5  L  +  G 


„  1  L5  ll’ 

Sum  =  +  —  1  “1 

4  fi5  4  G7 


L  3L°.5£t!i_L. 

‘V1'2  2G7  g5L+G 

5 

5  L  +  G 


‘ 3  74  P2 


+  1— 

G7 

3  l3  3  l[_ 
Sum  c  +  —  — r  ‘  7  7 

2  G5  1  G* 


-2  7  7 
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coefficient  of  cos  4g: 


9  2 _ 9^  2  :kl  -  2  .t! 

2°2  =  ’  64  g5  +  32  q7  64  Q9 


+ 


a2T2 


1  jJ  +  2  ll 


32  g7 


64  g9 


ii 

16  G 


7 

7 


Sum  =  0 


+  3^°2P2  =  0 

The  terms  with  cos  4g  disappear,  and  the  expression  for  F9*  becomes 

2  2 

after  substitution  of  the  appropriate  expressions  for  A  ,  AB,  B  and 
B  H2/G2  , 


2 

M  k2 

T  I  10 


3  L'6 

8  qi6 


15 

l2 

fl  -  2  h2 

-7HM 

32 

G'7 

V  G-2 

+ 


6.  2 
M  k2 

,,10 


1  -  16 


H' 


i,2 


+  15 


H 


„4, 


cos  2g’  , 


(29) 
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in  which  g  has  been  changed  into  g’,  which  is  permissible  in  a  part  that 
2 

has  kg  as  a  factor. 

5.  Secular  and  long -period  terms.  The  problem  is  now  reduced  to  the 
solution  of  the  system  of  canonical  equations  with  the  Hamiltonian 


F*  = 


A2 

L’3G’3 


♦  V- 


in  which 


F  *  -  F  *+F  * 

r2  r2s  +  *2p  > 


the  former  being  a  function  of  L',  G’,  H  only,  the  latter  depending  also 
on  g'. 

Let  S*  be  a  new  determining  function, 

S*  =  L"i'  +  G"g'  +  H"h'  +  Sj*  (L",  G”,  H",  g'); 


the  equation  F*  =  F**  may  then  be  written 

8S.  * 

Fn*  +  F,  *[  L",  G”  +  — i  ,  H" )  +F2s*  +  F2p*  =  F0**  +  F1**+F,** 


8g' 


This  equation  reduces  to 

8F.  *  8S,  * 

TT*-f**f*-it**  1  1  ,  t?  *  _  n 

*0  -  *0  ’  *1  -  *1  ’  — .  —  +  *2p  -  0  ’ 


9G"  8g’ 


(33) 


F  *  -  F  ** 

*2s  '  *  2 


(34) 
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Since 


8Y  3  M  k2  l\  .  5  H‘ 


8G"  2  l,3q„4 


•  Tf  2 


in  which  L'  has  been  written  for  L",  H  for  H",  equation  (33)  gives 


9g’ 


=  G" 


^2 


1/L-2  v*\rlmU jr+15  H 


®  VGm2  Qm4 


,„2 


4 

„4 


-1 


1-5 


h; 

Hm2 


cos  2g' 


Substituted  into 


G'  =  G”  + 


aSj* 

9g' 


there  results 

r 


G’  =  GM  { 


i  +  r2 


(\lL.v1 

Vg”2  g"4 


if i-n 

it*  2 . 


-1 


-5 


H 

ry*4 


1-5 


H‘ 

nti  2 


cos  2gM 
(35) 


Also 


sl*  =  G"y2 


L'2  _  L'4 

<  r  f  2  ^n4 


Ifl-llJL2^  -  I  h!  fl.5  h! 


-1 


16 


G"2/  2  g„4 


•m2 


sin  2gM 


and 


8S1*  L’ 

jr  =  r — -  =r'  +  y«- 

9L»  * 


if  i-n  H' 


-5 


H4  4.B5! 


-1 


.,,4 


•m2 


sin  2g", 


(36) 
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in  which  g’  has  been  changed  into  g". 

The  Hamiltonian  F**  is  a  function  of  L’,  G",  H  only.  It  is  given  by 
equations  (32)  and  (34), 

F**v  +  JLh-  (1 1  +  A  Ji"\  +  f2**  , 

2L’2  L,3G,,V  2  2  G "2/ 

where  the  last  term  is  the  first  part  of  the  right-hand  member  of  equation 
(29).  Let  nQ  be  defined  by 


Then: 
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Table  I.  Factors  arising  in  the  argument  of  the  perigee,  equation  (37) 


I 

25  cos^I 

4 

5  cos  I 

(1  -  5  cos2!)2 

1-5  cos2! 

0° 

1.562 

-1.250 

10 

1.539 

-1.222 

20 

1.476 

-1.142 

30 

1.395 

-1.023 

40 

1.350 

-0.890 

50 

1.552 

-0.801 

55 

2.140 

-0.839 

60 

6.250 

-1.202 

61 

10.58 

-1.576 

62 

25.72 

-2.381 

63 

234.76 

-6.956 

64 

115.72 

+4.716 

65 

12.45 

+1.120 

66 

3.790 

+0.792 

70 

0.232 

+0.165 

80 

0.001 

+0.005 

90 

0.000 

0.000 

This  completes  the  solution  of  the  problem  for  orbits  with  inclinations 
sufficiently  far  from  the  critical  inclination.  In  the  vicinity  of  the  critical 
inclination  the  results  may  become  illusory. 

The  data  given  in  Table  I  will  serve  to  show  that  for  an  orbit  with 
moderately  small  eccentricity,  no  significant  loss  of  accuracy  occurs  if 
the  inclination  is  even  as  close  as  1.58  from  the  critical  inclination,  63  26  . 
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The  second  column  In  this  table  gives  the  value  of  the  factor  that  am- 

tains  (1-5  cos2I)'2  as  it  appears  in  g',  multiplied  by  y2e  (1  -  e  )  •  1" 

h'  the  coefficient  is  -/Ttimes  that  in  g'.  No  term  of  this  type  occurs  in 

i\  The  third  column  in  the  table  indicates  the  magnitude  of  the  parts  of 

the  coefficients  that  contain  (1-5  cos2!)  1  as  a  factor.  In  X  and  g 

these  numbers  are  multiplied  by  ?  y2  times  a  factor  that  reduces  for 

e  =  o;  in  G'  and  h'  the  coefficients  are  diminished  by  a  function  of  e  that 
o 

contains  e  as  a  factor. 


6.  Comparison  with  results  obtained  by  Hill's  method.  For  e  =  o,  the 
motions  obtained  in  equations  (39),  (40),  (41)  may  be  compared  with  my 
results  obtained  by  Hill's  method  (Brouwer  1958).  With  L'/G"  =  1  these 
expressions  reduce  to 


diy 

dt 


=  n. 


dg” 

—  =n( 
dt  1 


dh”  n 
cos  I  —  =  n 

dt  ( 


*r2  (♦  £  - T  '°s¥*  Tf  C“‘r) 

y2  f-3cos2rj  +y22  ^+6  cos2r  -  cos4r^ 
l+y2  f-S+ncos2!"^  +y22  |+ ~ -36cos2r  +  2£2cos4rj 


d(i"+g") 
dt  c 


The  sidereal  mean  motion  is 

d(l"+g")  tm  dh” 

n  =  v  .  ?J.  +  cos  I  — 

s  dt  dt 


(42) 


ns  =  no 


l  +  y2  I-3  +  9COS2!"'  +y22  -30cos2r  +  ^.  cos4!’ 


or 


no  =  ns  1+y2  f+3'9cos2r'j  +y22  (+7  -  24cos2r  +  ^cos4rj 


Thus 


«N?.S  1  +y2  ^+3cos2rj  +y22  ^3cos2r+jcos4r) 

CosI"^"=nB  y2  ^ScosVj  +  y22  3  cosV  - 1  cos4^  .(43) 

In  the  Hill  method  applied  to  the  artificial  satellite  problem  I  found 
(Brouwer  1958,  p.  435) 

d(iM+g,t)  =  1  +  3y*2cos2lH  +y2*^(~2  "  30cos  9  ^ 

dt  \  “  ' _ 

n  „  and  I„  being  used  to  distinguish  these  symbols  from  ng,  I  in  the 
sH  H 
present  article. 

The  expression  in  square  brackets  is  equivalent  to  the  value  of 
1/(1  +  p)  for  y4  =  o  with  y2  replaced  by  y2*.  In  order  to  compare  with 
the  result  obtained  above  it  should  be  noted  that  y2*,  cos  IH  and  ngH  differ 

from  y2,  H/G"  and  ng. 

In  the  present  theory,  if  L'  =  (paQ)  '  , 


p2^ 2  ^2  n  _  p2  -  P2 

~4  2  n°  .,3 

L  a0  2 

a0 


In  the  Hill  theory 


’  nsH 


Hence 
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v  *  2 
y2  a0 


=  1+4 y2  (1-3  cos^)  , 


from  equations  (42)  and  (44).  Hence 
yg*  =  ?2  '  ^2  ^  "  3  cos2Ih)  . 

The  inclination  1^  in  the  Hill  method  was  defined  as  the  maximum 
latitude.  In  the  present  theory  this  does  not  apply  to  the  inclination  ob¬ 
tained  from  cos  I”  =  H/G”  in  view  of  the  presence  of  the  term  of  the  first 
order  in  y ^  and  zero  order  in  e  in  G/G'  For  e  =  o  the  relation  is 


G  =  G’  1  +  -  y9  sin2 1  cos  (2g  +  2  i)  . 
2 


Then 


A(coszI)  = 


2 

_  2H_  AG  =  _3  sin2Icos2Icos(2g+  21) 

2  G’  L 

G’^  U 


in2 1  =  sin2I"  1+  3y,  cos2!"  cos  (2g  +  2 1)  . 


sin2/3  =  sin2Isin2(g+A)  =  sin2I"  1 +3y2cos2I"cos(2g  +  2A)  sin2(g+l) 


This  is  a  maximum  for  2g  +  2 1  =  180“,  when 


sin20  =  sin2  IH  =  sin 
Hmax  H 


in2r'  l-3y2cos2r  , 
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from  which 


COS^Ijj  =  cos2!" 


1  +  3y2  sin2 1" 


(46) 


The  introduction  of  (45)  and  (46)  into  (44)  yields 


dt 


sH 


l  +  3y2cos2r  +  y22  (-  cos2 1"- 3 cos4 r 


.(47) 


Finally 


dh” 


n  -  d(l"  +  g")  +  ensf  &  ,  n  -  d(1"±^  +  cosIH  -  , 

sn  dt  n  dt 


dt 


dt 


or 


n_„  -  no  =  (cos  IH  -  cos  r )  -^’  =  -?  y2  sin2 1"  cos  I"  — 


sH 


dt 


=  -  1  n  v22sin2  I”  cos2r  . 
2  s  L 


Introduction  into  (47)  yields 


d(lM  +  gM)  =  n 


dt 


1  +  3y2  cos2r'  +  y22  (s  cos2r’  +  -2  cos4!’^ 


in  agreement  with  the  expression  for  obtained  previously. 

In  addition  to  establishing  the  agreement  between  the  results  for  e  =  o 
obtained  by  two  different  methods,  the  discussion  serves  as  an  illustration 
of  the  caution  necessary  in  the  comparison  of  results  to  the  second  power 
of  y2  obtained  by  different  treatments  of  the  same  problem. 
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7.  The  fourth  harmonic.  The  coefficients  of  these  terms  are  so  small 
that  it  is  unnecessary  to  consider  the  terms  of  short  period.  Hence  the 
contributions  to  the  Hamiltonian  m^y  be  limited  to  the  parts  independent 
of  J,  and  only  the  first  powers  are  needed. 

The  expression  for*U^  is 

U,=-4fl-10  sinV-^sinV). 


Since 

sin2/3  =  —  sin2I  [  1  -  cos  (2g  +  2f)l  , 

2 

sin^jS  =  —  sin^I  [3-4  cos  (2g  +  2f)  +  cos  (4g  +  4f)l  , 
8 

this  becomes 

g  - 

=  —  /l  -  li  cos2!  +  M  eosV)  al 


10  V8 


(  5  20  2t  35  4A  a5 

+ - +  —  cos  i  -  — cos  I  —  cos  (2g  +  2f) 


f 35  35  1.  35  4?\  a5 

+  —  -  —  cosT  +  —  cos  I )  — cos  (4g  +  4f)  . 
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The  secular  part  is 


3_  -  15^  +  35  5_  3 


l3g7  IV 8  4  g2  8  g4 J  \2  2  a 


+  ?0  h!.35  HlVl.  3  GiVos%'  _ 

V.  6  3  g2  6  Q4y  \4  4  tf)  J 
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The  part  independent  of  g  is  added  to  ^g*  > 


*4F2s*  - 


n\ 


f 15  L^_ 

V16  G'7 


_9 

16g.5 


1  -  10 


> 


the  part  having  cos  2g  as  a  factor  is  added  to  F2p 


*  5  *^4 

A4F2p  &  ~To 


P_ 

m  7  (V  ^ 


8^+7  — ^  cos  2g\ 
G'2  G'V 


Now 


9 

9g 


2  L,3Gm4 
2  n\ 


4V1  v* 

G"  / 


By  integration, 


K  2  k 

A/S,  *=  — -H - —  G' 

4  1  24  L, 4  k2 


If  now 


y4  = 


> 


1 

cos  2g’ 


cos2g’  . 


sin2g?  . 
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so  that 

A 4  =>4 
h\  r2 


it  follows  that  the  additions  to  the  right-hand  members  of  equations  (35) 
to  (38)  are 
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The  additions  to  the  secular  motions  are 


8.  The  third  and  fifth  harmonics.  Odd  harmonics  are  to  be  considered 
only  if  symmetry  of  the  earth  with  respect  to  the  equator  does  not  exist. 
A  study  of  the  motion  of  1958  02  led  to  the  introduction  of  these  odd  har¬ 
monics  by  O’Keefe,  Eckels  and  Squires  (1959). 

The  expression  forUg  is 


u 


pA 


3.0 


'5  .  30  3  •  n 

sin  /3  -  —  sin 


With 


sin**/} 


sin3 1 


—  sin  (g  +  f )  -  —  sin  (3g  +  3f) 
4  4 


this  becomes 


TJ,  = 


ma 


3.0 


1  sin  1+  —  sin0 1 )  sin  (g + f)  -  1  sin3 1  sin  (3g  +  3f) 

2  8 
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£  J_ 

2  2  ^  2y  2v2 

For  the  sake  of  brevity,  let  e’  sin  I’  stand  for  (1  -  G’  /V  )  (1  -H  /G’  /  , 
eM  sin  I”  for  the  same  expression  with  G’  replaced  by  GM.  The  secular 
part  may  then  be  written 


A3F2p* 


3,0 

L'3G'5 


e’  sin  I’ 


l5  sin  g'  . 

8  g'2/ 


There  is  no  part  independent  of  g\ 

It  follows  that 

A  2 

—  AqS«*=  —  — —  — _  G”  e”  sin  I"  sin  g’, 

8g’  4  ^,2  kg  q.«2 


1  ^  Q  A  t  1 2 

*  ^  G"  e"  sin  I"  cos  g1  . 


A  =  +  -  -ii - 

4  L’^  k2  G"2 


If 


r3 


^A3.0 


a3g’  = 


1  h _  G1  eM  sin  I"  sin  g”  . 

4  y2  G”2 


With 


JL  (e  sin  I) 
3L 


sin  I  G2 


J!_  (e  sin  I)  = 
8G 


sin  I  G  +  e  H2 
e  ^2  sin  I  q3 
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JL  (e  sin  I)  =  -  — —  —  , 


8H 


sin  I  q2 


the  following  results  are  easily  obtained: 


v  =  cosg„ 

6  4  y2  e"  L' 


4  r2  G"2  V  e"  sin  r  G"2. 


1  r3  L’2  H  e" 


Ajh1  =  +  ±  _  _ 


4  y2  G"2  G"sinI" 


cos  g"  . 


The  fifth  harmonic  is 


T7  *^5.0 
u  5  =  - 


6  L 


15  sin  0  -15  sin3  j3  +  l3  sin5  0 

8  4  8 


With 


sin5 


sin5 1 


1  sin  (g  +  f)  -  —  sin  (3g  +  3f)  +  —  sin  (5g  +  5f) 
8  16  16 


this  becomes 


u*  = 


5.0 


15  sin  I  -  —  sin3!  +  —  sin5l)  —  sin  (g  +  f) 


16 


64 


6  6 
+  sin3 1  -  —  sin5 1^— sin  (3g  +  3f)  +  —  sin5 1— sin  (5g+  5f) 


16 


128 


128 
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w 


M  As-0  G"  e”  sin  I” 
L,8k2  G"6 


With  the  aid  of 


*5 


^^A5.0 
L' 10 


the  resulting  long -period  terms  become: 
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I 


-fca«£j 


mmtmammmm 


37 

9#  Formulas  for  computation.  For  convenience  of  computation  the  per¬ 
turbations  in  the  Keplerian  elements  a,  e,  I  are  given  instead  of  those  in 
L,  G,  H. 

The  adopted  force  function  is 


1  - 


3  sin2/3 


1  -  10  sin  j3 


1  sin/3  +1  8in3^J 


^fesW-™sinb^8in*{i 
8  4  8 


in  which  k£  is  a  small  quantity,  and  k^,  A^  q,  Ag  q  are  assumed  to  be  of 
order 

The  secular  motions  have  been  computed  to  Ofcg2),  the  coefficients 
of  periodic  terms  to  Ofkg). 

Basic  constants: 

a"  =  semi -major  axis  constant 
e”  =  eccentricity  constant 
I"  =  inclination  constant 

n  =  (i1/2  a"‘3//2  17.04337  (a"/R)'3'/2  rev./day 

R  =  equatorial  radius 
Abbreviations: 
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r4 
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-10 

r5v 
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It  is  customary  to  use  for  the  second  harmonic  the  coefficient  J; 
Jeffreys  (1954)  used  for  the  fourth  harmonic  the  coefficient  D.  The  re¬ 
lations  between  J,  D  and  y2,  >4  are 


3  n/R' 

r4  =  i?Dt, 


Strictly  speaking,  e"  +  6^,  6’  =  cos  (r'  +  djD.n’  =[l  -  (e"  +  6le)‘ 
should  be  used  in  the  computation  of  the  periodic  terms,  but  since  the 
short-period  terms  are  obtained  to  0  (4),  it  is  of  no  consequence  if  con¬ 
tributions  of  O  (k2)  are  omitted  in  expressions  that  have  y2  as  a  factor. 
Similarly,  A" ,  g"  might  be  used  in  computing  f,  r';  but  since  i  ,  g  a.  e 
available,  their  use  does  not  complicate  the  calcula  ion. 

The  formulas  are  applicable  for  any  eccentricity  e  1  and  any  in¬ 
clination  with  the  exception  of  inclinations  near  the  critical  inclination, 

for  which  1-5  cos2I  appears  as  a  small  divisor. 

The  appearance  of  e"  as  a  divisor  in  the  short-period  terms  in  e  is 
apparent  only.  The  expressions  that  are  multiplied  by  e"*  contain  e"  as 

a  factor,  either  implicitly  or  explicitly. 

In  the  short-period  terms  in  Jl  and  g  a  divisor  e"  occurs  also,  but 
for  the  calculation  of  the  position  only  g  +  I  +  equation  of  the  center  is 
needed.  In  g  +  X  the  divisor  e"  is  not  present. 

Singularities  in  some  of  the  elements  also  occur  for  very  small  in¬ 
clinations;  again,  no  singularity  is  present  in  the  coordinates.  In  such 
cases  it  may  be  found  convenient  to  modify  the  formulas  and  obtain  ex¬ 
pressions  for  the  perturbations  in  coordinates. 

Secular  terms 

1"  =  ’’mean"  mean  anomaly 

= v- 1  +  \  y'2v  1  + 


1/2 


3  ,  2„ 

+  —  y9  v 

32 


2\„4 


- 1 5  +  16tj  +  25i}2  +  (30  -  96q  -  90)]2)  6  2  +  (105  +  144ij  +  25 7)  )6 


15  ,  q„2 

+  16  y4  ^ 


3  -  306 2  +  3 56 4 


+  V 
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g”  =  mean  argument  of  perigee 

■  V  jj  A  <-  1  *  m2> 


3  ,  2 

+  alr2 


-  35  +  24 11  +  25jj2  +  (90  -  192j?  - 1 267J2)  9  2 + (385  +  360?) + 45I}2)  fl‘ 


+  -y,4  1^1  -  9»?2+ (-  27°  +126tj2)02+(385-1891)2)8^]^  +g0" 


h"  =  mean  longitude  of  ascending  node 
r 

=  "o1  <  -  3>,2  9 

+  3  /  2  ["(-  5  +  12r)  +  9J)2)  9  +  (-  35  -  36t)  -5t;2)  93J 
8  L  L 

+  -  r'4  (5  -  8 1)2)  9  (3  -  W2)  y+  h0" 


Long-period  terms: 
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h  =h’  -  -Iy'0X 
2  L 

x  ^6  (f '  -lf+  e"  sinf)  -3sin(2g’  +  2f)  -3e”  sin(2g’  +  f)  -e”  sin  (2gf  +  3f’)J  . 
f’,  r’  are  to  be  computed  from 


1 

2 

iL  sinf  =(l-e"2)  sinE' 
a" 

__  cos  f  =  cos  E’  -  e” 
aM 

=  1  -  e"  cos  E' 
a" 


For  the  calculation  of  the  coordinates  at  any  time  the  complete  values 
of  e  and  X  should  be  used  for  the  solution  of  Kepler’s  equation, 

E  -  e  sin  E  =  X 

and  subsequently  r  and  f,  which  may  then  be  used  in  the  formulas: 
x  =  r  [  cos  (g  +  f)  cos  h  -  sin  (g  +  f)  sin  h  cos  I] 
y  =  r  [  cos  (g  +  f)  sin  h  +  sin  (g  +  f)  cos  h  cos  I] 
z  =  r  sin  (g  +  f)  sin  I 
A  convenient  alternative  form  is: 
x  =  A„  (cos  E  -  e)  +  B  sin  E 

X  X 

y  =  A  (cos  E  -  e)  +  B  sin  E 

J  J 

z  =  A  (cos  E  -  e)  +  B  sin  E 
z  z 


Ax  =  a  [  cos  g  cos  h  -  sin  g  sin  h  cos  I] 
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B„.-. 


It-/) 


sin  g  cos  h  +  cos  g  sin  h  cos  I] 


A  =  a  [  sin  g  cos  h  cos  I  +  cos  g  sin  h] 


1 

2 

[  cos  g  cos  h  cos  I  -  sin  g  sin  h] 


Az  =  a  sin  g  sin  I 


B. 


=  a 


1_ 

2 

cos  g  sin  I 
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